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INTRODUCTION 

This final report is an overview of the work done on Cooperative Agreement 
NCC 3-55 with the Solid State Technology Branch of the NASA-Lewis Research Center 
(LeRC). Over the period of time that the agreement was in effect, the principal 
investigator - and in the last three year - the co-principal investigator (Walter*Rr 
.Lambneefot) worked on a significant number of projects and interacted with members 
of the Solid State Technology (SST) branch in a number of different ways. For the 
purpose of this report, these efforts will be divided into five categories: 1) work 

directly with experimental electrical transport studies conducted by members of the 
SST branch; 2) theoretical work on electrical transport in compound semiconductors; 
3) electronic structure calculations which are relevant to the electrical transport in 
polytypes of SiC and SiC-AIN alloys; 4) the electronic structure calculations of polar 
interfaces; and 5) consultative and supportive activities related to experiments and 
other studies carried out by SST branch members. Work in these categories will be 
briefly discussed b’feTow: Much of the past work was discussed in more detail in the 
various proposals for this grant. 

Copies of references [4], [5], [6] and [7] are attached to this report. Copies 
of references [1], [2] and [3], which also acknowledge support by this grant, are not 
included as the monitor of this grant is a co-author of those publications and thus has 
copies of them. 
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WORK CARRIED OUT UNDER THE AGREEMENT 


(1). Work Directly Connected with Experimental Electrical 
Transport Studies Conducted by Members of the SST branch: 

Because of SiC potentialities in a wide range of applications, in particular, 

t 

semiconductor applications at high temperature, high fields and high frequencies, this 
material has attracted considerable interest over the last decade. In the early 1980s, 
researchers at LeRC were the first to grow epitaxially large area crystal films. This 
was done by the chemical vapor deposition (CVD) method. Thus, there was 
considerable interest at LeRC to determine the electrical transport properties of 
those films. We note here that while SiC can form in a large number of 
crystallographic structures (polytypes), those films had the 3C (the cubic) structure. 
The electrical properties are obtained by so-called Hall measurements. A series of 
such measurements on several films produced at LeRC was initiated by Drs. E. J. 
Haugland and S. A. Alterovitz of the SST branch. The PI worked with these 
researchers to interpret the data in terms of the band structure of 3C SiC. The 
results of these investigations appeared in Applied Physics Letters [1]. 

About the time of this work, two other groups of researchers published the 
results of measurements on similar SiC films. Their analysis lead to different 
interpretations than ours. Our analysis of their results, which was also published [2], 
backed our interpretation, which is now the widely accepted interpretation. 
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Another project in this category was also related to experimental work by Drs. 
Haugland and Alterovitz along with that of a visiting researcher, Dr. P. P. Szydlik. This 
involved the study of the properties of the two-dimensional electron gas in a ln 0 .iGa 0 .9 
As/GaAs heterostructure. In this structure the electrons reside in the InGaAs alloy, 
which has a smaller band gap than the GaAs. The experiments consisted of Hall and 
Shubnikov de Haas (SdH) measurements. The latter involves oscillations in the 
magnetoresistance. From analyses of the SdH oscillations as a function of 
temperature and as a function of magnetic field, it is possible to determine an 
effective mass m* and an electron scattering time x SdH- From the Hall mobility - 
determined in the Hall measurements - another electron scattering x H could be 
determined. It was found that the ratio x H /xsdH was about 5 to 6. The fact that the 
ratio was considerably larger than one was explained by the fact that while small 
angle scattering does not contribute appreciably to x H , it does contribute to xsdH- 
Relevant scattering mechanisms which could contribute to the difference were 
identified and the value of m* was discussed in terms of known values for the pure 
compounds. A paper on this work was published [3]. 

(2). Theoretical Work on Transport in Compound Semiconductors: 

The aim of this work was to set up a reliable means for calculating the 
transport properties, in particular, the Hall mobility and the so-called Hall coefficient 
r H (T) (which is given by Rh = r H /en where R H is the Hall constant and n(T) is electron 
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density). Both quantities, especially (i H . depend on the various scattering mechanisms 
that are operative and on the nature of the electronic energies in the vicinity of the 
relevant band extremum, e.g., conduction band minimum (or minima). For SiC, and 
other semiconductors of current interest, both the mechanisms and the details of the 
band edges are not known with accuracy and certainty. The goal is to compare the 
calculated m-h(T) with suitable experimental results over a reasonable range of T and 
thereby pin down the relevant scattering mechanisms and certain aspects of the band 
edges. 

The PI, with the help of two CWRU undergraduate physics students, has set up 
programs to numerically solve the Boltzmann transport equations. All the principal 
scattering mechanisms relevant for compound semiconductors are included. The 
programs have been checked extensively. Of most immediate interest to us is the 
high temperature mobility of 3C SiC. Subsequent to that will be the study of the 
mobilities in the other polytypes; that will probably be a somewhat more involved 
undertaking. 

(3). Electronic Structure Calculations Relevant to Electronic 
Transport in Polytypes of SiC and Related Materials: 

A knowledge of the band structure is important for the transport because 
transport is effected by the nature of the band edge: whether it is degenerate or 
nondegenerate, whether it is at the center of the Brilloun Zone (BZ) and hence is a 
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single extremum, or is at some other k value in the BZ and hence is a “multi-vallied” 
extremum. In addition to having a large density of states, a multi-vallied band edge 
introduces the possibility of “intervalley scattering”. Also important is the particular 
location of the extremum and the effective mass of the band edge. 

The study of the band structure(s) of SiC is a significant project because, as 
noted earlier, SiC forms in numerous different structures, or polytypes. The different 
structures can be described by different orderings of the close-packed layers of 
atoms, these layers being the (111) plane of the cubic or the (0001) planes of the 
various hexagonal or rhombohedral structures. For our present purposes the 3C, 2H, 
4H and 6H polytypes are the most important. The latter two are the most common 
and easily grown forms of SiC while the 2H (wurtzite) form is difficult to grow. The 
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various structures differ in the frequency and ordering of the “cubic” ( JF ) and 
“hexagonal” (» stackings of the successive layers. With 3C having all the cubic and 
2H all the hexagonal stacking, the sequence of the four polytypes in terms of 
increasing “hexagonality” is 3C, 6H, 4H and 2H. 

Since the position of the valence band maximum of all the binary compound 
semiconductors, with the exception of AIN, are at k=0, it is the conduction band 
minima which are of most concern. (All the values of the parameters describing the 
degenerate valence band maxima are unknown and thus, are of interest.) We first 
studied the 3C and 2H, the simplest forms of SiC. In all of our self-consistent 


calculations, we obtain the “energetic” properties (lattice constants, bulk moduli, and 
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cohesive energies) in addition to the band structure. It was found that the band gaps 
were indirect in both structures with the minimum for 3C being at the point X (of the 
fee BZ) and that for 2H being at point K of the hexagonal BZ with a state at M being 
close in energy. We note that the points K and M are not the hexagonal equivalents of 
the point X when 2 BZs are appropriately superimposed. In addition to the “shifting” 
of the position of the conduction band edge, the band gap for 2H was found to be 
almost 1 eV larger than for 3C. To obtain a qualitative understanding of the above 
results, calculations for the closely related compound semiconductors AIN and GaN 
were carried out. An understanding was achieved and was described in a publication 
[4] acknowledging this grant. 

More recent work included calculations of the 4H and 6H polytypes. With the 
results of these studies, we were able to determine the trends in the various band 
gaps including the minimum gaps of the aforementioned four polytypes as a function 
of hexagonality. This among other things was addressed in a publication [5] 
acknowledging the grant. 

Further work studying the details of the conduction band edges is in progress. 
We note that recent transport studies conducted at the CREE Corporation have 
shown that the mobility of 4H SiC is higher and much more isotropic than that for the 
6H form. These results can have significant implications for device applications, and 
represent one “hottest” unsettled issue in the SiC community. Tentative results from 
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our calculations may have some bearing on this question and will be pursued 
vigorously in the near future. 

Finally, we note the work on alloys of SiC and the closely related compound AIN 
which is closely latticed matched to SiC. These alloys are of interest because AIN has 
a significantly larger band gap than SiC (for wurtzite AIN, E g = 6.3 eV) so that the 
alloys will allow us to “tailor” the band gap over a large range. Moreover, for the 
equilibrium (wurtzite) form, AIN is direct. Thus, there is a crossover from an indirect 
to a direct band gap. A direct gap is advantageous for optical emission. Work on 
these alloys was briefly discussed in a report to the 1994 International Conference on 
SiC [6]. A more complete report will be written shortly. 

(4). Electronic Structure Calculations of Polar Interfaces 

Since semiconductor heterojunctions are now playing an important role in 
electronic technology and in basic physics, the study of the electronic structure at 
interfaces between semiconductors hardly needs justification. Such studies are 
generally done using the supercell approach, i.e., by studying a periodic superlattice of 
the two constituent semiconductors with a suitable (generally large) unit cell. 
However, when the interfaces are charged (the condition which characterizes “polar” 
interfaces), the calculated interfaciai energy is the sum of the energies of the two 
oppositely charged interfaces. Since two interfaces are inequivalent, it is not possible 
to obtain the energies associated with the individual interfaces. Similarly, because of 
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the existence of an electrical field associated with the interfacial charges, it is also not 
possible to separate out the band offsets at the two interfaces. A scheme for 
circumventing this problem was worked out and published [7], The approach was 
applied to polar inversion domain boundaries in GaAs and NaCI. 

(5). Consultative and Supportive Activities 

The PI has consulted with several members of the SST branch in connection 
with some of their projects. We briefly discusses a number of interactions. 

The most frequent interactions were those with Dr. E. J. Haugland. The PI had 
numerous discussions with him about anticipated transport studies particularly those 
for SiC at relatively high temperatures. There were also numerous discussions with 
Dr. Haugland about his recent Hall data on 4H and 6H films. 

The PI has, over the past several years, sought to obtain samples of 3C SiC of 
reasonably good quality for Dr. Haugland to study in the high temperature Hall 
apparatus that he has set up. Recently, we have finally succeeded in obtaining such 
samples from a research group in the Naval Research Laboratory. Hopefully, this will 
allow a definitive determination of the high temperature mobility of the cubic form. 
This will be compared with our theoretical calculations. 

There were also discussions with Dr. Haugland about the interpretation of his 
electrical transport measurements of a number of GaAs crystals grown in space (on 
the STS-40 flight, June 1994). The properties of these samples differed from those of 
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earth-grown samples, which were otherwise as similar as possible to their space- 
grown counterparts. Dr. Haugland has reported on the data and on the (tentative) 
interpretation of them in internal NASA meetings 

The PI also consulted with Rafael Mena on the changes in the electrical 
transport properties of the two-dimensional electron gas at the interface of 
AIGaAs/InGaAs heterojunctions produced by illumination. Among other issues, the 
effect of electron occupancy of the second and possibly higher subbands was 
discussed. 

Finally, there were several discussions with Dr. David Rosenfeld on some 
transport problems in SiGe alloys in certain transistors. 
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ELECTRONIC STRUCTURE AND RELATED PROPERTIES 
OF TETRAHEDRALLY BONDED WIDE-BAND-GAP MATE- 
RIALS CONTAINING EARLY ELEMENTS OF THE PERIO- 
DIC TABLE 


W. R. L. LAMBRECHT. C. H. LEE. Iv. KIM. A. G. PETUKHOV. 

E. A. ALBANESI AND B. SEGALL 

Department of Physics 

Case Western Reserve University 

Cleveland. OH-44106-7079 


1. Introduction 

The tetraliedrally bonded materials involving early elements of the periodic 
table, such as Be. B. C, and N have rather extreme properties compared 
to the conventional tetraliedrally bonded semiconductors. These unique 
properties make them interesting for a variety of applications, including 
their use as hard coatings for mechanical tools, metal/ ceramic compos- 
ites. heat-sinks, and electronics. The latter are the most demanding type of 
applications because they require semiconductor grade purity single crys- 
talline material in order to realize the superior performance these materials 
promise to have. But their extreme properties also makes them difficult 
materials to grow in the form of pure single crystals. Their properties are 
thus not yet very well known. In order to make a rational choice of material 
for specific applications, it is necessary to understand the interrelationships 
between the intrinsic materials properties and the trends in these properties 
with atomic number, crystal structure, and so on. Virtually all materials 
properties can in ultimate instance be related to the underlying electronic 
structure. The purpose of this paper is to describe the basic trends in elec- 
tronic structure of these materials and some of their related properties. 
The results presented here are based on a systematic study of the elec- 
tronic structure of these materials carried out over a number of years in 
our research group. Although many of the results discussed here were pre- 
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seated elsewhere they were 110 c included in an overall description ol the 
trends as presented here. Some new results are included as well. Extensive 
references are provided as a guide to our previous work not covered here in 
detail. 

We first describe some of the relationships between intrinsic properties 
and their importance for various types of applications. Next, we briefly de- 
scribe our computational approach. Then, we present calculated results for 
the groundstate properties such as the equilibrium lattice constants, bond 
strength, elastic constants, as weil as the underlying electronic structure — 
in particular the band gaps — for a number of materials systems. 

We classify the materials as follows: 

- materials containing only elements from the second row of the periodic 
table: diamond. BN. BeO 

- related Be compounds: BeCN? and Be 2 C 

- poly types of SiC: 

- B-compounds: BN. BP. BAs 

- the group Ill-Nitrides: BN. AIN. GaN. InN: 

- related II-IV-N 2 materials: 

- alloys among III-N’s and heterovalent alloys. 

2- Relations among intrinsic properties 

Wide-band-gap semiconductors have a number of advantages for electronic 
applications over conventional semiconductors. Among these are larger 
breakdown fields, resistance to harsh environments, ability to operate at 
high-temperatures, and high-thermal conductivity. These various proper- 
ties are important for a variety of applications where conventional semi- 
conductors fail. For example, in conventional semiconductors, intrinsic car- 
rier concentration at a few 100°C wipe out the effects of doping and make 
semiconductor devices inoperable. Obviously, since intrinsic carrier gener- 
ation varies exponentially with the gap, wide band-gap materials have an 
important advantage in this respect. Intrinsic carrier generation also plays 
an important role in the time-decay of space charges induced by optical 
excitation of e.g. an MIS or Schottky diode. This may enable new kinds of 
semipermanent memory and UV detetectors. 

At high temperatures, the formation of various types of defects also 
start to degrade the material. For example, dislocation movement proceeds 
by kink formation and migration. Vacancy and interstitial production also 
have a thermal activation energy. All these processes require bond breaking 
and thus the stronger the covalent bond the more resistant the materials 
are to high-temperature degradation of the crystalline quality. The wide- 
band-gap materials we will consider here derive their wide gaps from strong 
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covalent bonding and not exclusively from hiah ionicity. They are thus 
also strongly bonded. This is also what gives these materials resistance to 
radiation damage. 

The above effects are not only important for higli-temperature operation 
but also for high-power and high-frequency applications. In fact, both oi 
these lead to generation of heat. The thermal conductivity is thus another 
important performance parameter for a semiconductor. The wide-band- 
gap tetrahedrally bonded materials have a high phonon-mediated thermal 
conductivity because of their high sound velocities, which in turn are closely 
related to the stiffness of the bonds. The latter can best be described in 
terms of bond strengths, elastic constants and the force constants involved 
in certain phonon modes. 

Another important parameter for high-speed electronics is the static 
dielectric constant. This is because in a very general sense, electronic devices 
are based on the transfer of charges from one spatial area to another and 
thus involve capacitance. One clearly wants as low a dielectric constant as 
possible. Also good insulators with low dielectric constant are required for 
isolating various active parts of a device from each other. Thus, as one starts 
to use wider gap semiconductors, one also needs even wider gap insulators 
for packaging applications. The electronic dielectric constant at frequencies 
above vibrational frequencies is roughly inversally proportional to the 
band gap. In partially ionic materials, the static dielectric constant involves 
also a contribution from the ionic motions. Thus ionicity is an important 
parameter. Ionicity can be defined in terms of the ratio of the atomic energy 
level difference to the effective gap between bonding and antibonding states. 
Phillips' ionicity /, [1] is related to Harrison's [2] polarity ap by /,■ = a z P 
with 

a P = | E c - EJ\/J\E C - E a \ 2 + 4ti 2 . (1) 

Here E c and E a are cation and anion atomic energy levels and i3 is the 
covalent interacti on between sp 3 o rbitals pointing towards each other. The 
quantity Eq — \J\E C — E a p -I- 4/J 2 can be interpreted as an average effective 
gap between cation and anion levels. Effective values for these levels can 
be extracted from our first- principles calculations [3] and define an ionicity 
scale. The ionicities of some of the materials discussed here are given in 
Table 1. 

The underlying reason for the unique properties of the materials involv- 
ing early elements of the periodic table can be understood on the basis 
of the simple molecular bond picture embodied in Eq.(l). Because of the 
absence of lower energy p-like orbitals in the core, the 2p levels are anoma- 
lously deep compared to 3p levels. Also. 2s levels feel the nucleus more 
strongly than 3s electrons because the nucleus is less effectively screened 
by the small core consisting only of Is electrons. This relatively compact 
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TABLE i. Polarities calculated from 

LMTO-ASA calculations for equal cation and an- 
ion sphere sizes 


SiC 

0.473 

BN 

0.475 

AIN 

0.807 



BP 

U.018 

GaN 

0.771 



BAs 

0.019 

InN 

0.792 

BeO 

0.818 

BeC* 

0.588 

BeN* 

0.701 


Individual bond polarities in BeCN?. 


nature of the 2s and 2p orbitals is responsible for the small lattice constant 
and strong covalent bonding. The strong covalent interaction d leads to a 
large effective gap Erj. 

When atoms of the second row are combined with atoms of the third 
row. there is also an anomalous difference in atomic energy levels. This 
leads to a strong ionic component to the bonding for SiC (even though the 
latter is a IV- IV compound) and AIN and a surprisingly low ionicity for 
BP. The latter occurs because the 2p level of B is sufficiently deep that 
its position approaches that of the P 3p level. By Eq.(l) that implies low 
ionicity. In fact, taking into account that the effective charge on the V- 
element is given by q = 1 - 4ap, one finds that it-actually has positive 
charge when ap < 1/4. This is the case for BP and BAs which as a result 
have inverted ion character: P or As behave as cation while B behaves as 
anion despite the fact that normally the group-III element behaves as cation 
and the group-V element as anion. On the other hand, the enhancement 
of ionicity due to the deep 2s and 2p levels for SiC and III-N compounds 
contributes significantly to their large gaps. 

The relative extent of 2p versus 2s orbitals compared to those of 3s 
and 3p has another important effect. It implies that there is less advan- 
tage in promoting a 2s electron to a 2p orbital to form directed bonds. In 
fact, strong bonds can also be realized by means of 7r bonding between 2p- 
orbitais directly. This means that for elements from the second row of the 
periodic table there tends to be a competition between tetrahedral bonding 
and planar sp 2 plus n bonding. Thus, C and BN have competing layered 
structures: graphite and h-BN. BeO does not exhibit such structures [4] 
which is related to the higher ionicity of tills II- VI compound. Indeed ion- 
icity tends to favor higher coordination. 

Ionicity also plays a role in the relative stability of cubic versus hexag- 
onal stacking of tetrahedraliy bonded layers [5], Thus, the relatively low 
ionicity compounds BP. BAs and BN and (zero ionicity ) elemental C' fa- 
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• or the cubic structure while the Ill-nitrides have wurtzitic structures and 
SiC exhibits polytvpism. lonicitv plays also an important roie in electronic 
transport. The polar mode scattering increases with increasing effective 
charges. 

Wide band gaps have another important application in optoelectronics 
where they extend the range of light emission frequencies into the blue and 
UV. In this context, it is important whether the gap is direct or indirect. 
Other important parameters for electronic applications are effective masses, 
deformation potentials, and other band structure details. Some examples 
will be discussed below. 

Finally, an important aspect of potential semiconductor materials is the 
possibility of doping the material n and p-type. Diamond has proved to be 
quite difficult to dope n-tvpe. GaN was difficult to dope p-tvpe although 
effective p-tvpe doping with Mg was recently achieved using post-doping 
annealing and Low energy electron beam irradiation iLEEBIl treatments 
[6]. AIN has sofar resisted effective doping of either type. BN has been 
reported to be both n and p- type dopable [7]. Although this is an important 
issue, we will not discuss it further because we have not done studies of this 
aspect in our group. 

In summary, we have here given some examples of how the electronic, 
structural and bonding properties are interrelated and how they impact 
device applications. In the remainder of the paper we focus on the electronic 
structure. 

3. Computational Method 

The theoretical framework of our first-principles calculations is density 
functional theory [8]. This approach basically reduces the many- body prob- 
lem of interacting electrons to that of independent electrons in an effective 
potential with the same ground state electron density as the real system. 
The total energy of the system is expressed in terms of the Couiombic en- 
ergy of the electron charge density and the positively charged nuclei, the 
kinetic energy of the independent model electron system and a remainder 
which is called the exhange-correiation energy. The method is in principle 
restricted to the study of groundstate properties. The local density approxi- 
_ mation (LDA) allows one to incorporate effects of exchange and correlation 
in a mean-field approximation. It assumes the same functional form for the 
exchange-correlation energy as in the homogeneous electron gas but ap- 
plies it locally to each volume element of the real inhomogeneous system: 
Etc = / e rc [n(r)]n(r)d 3 r. We use the parameterization of the exchange and 
correlation functionals by Hedin and Lundqvuist [9]. 

An important point to realize about this theory is that the eigenvalues 
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of the effective one-eiectrou Schrodinger equation 


- — A 4- i’tf(r) + t , re lrjtt , n fc(r) = ( n kVnk{ r). 
2 m 


where r// is the Hartree potential combining the electron and nuclear eiec- 
trostatic potential, and v rc is the exchange-correlation potential, do not 
have the meaning of single particle excitations. For a periodic solid the 
( n k constitute the band structure. In other words the band structures we 
calculate in this theory form a useful construct to understand total energy 
properties and bonding but do not represent the energies required to ex- 
tract electrons from the system as measured for example by photoemission 
or as involved in optical excitations (excluding excitonic effects). The true 
quasiparticle excitations are given be the related equation: 


(“T^-A + i'//(r)] #„*(!•) -F j Z:. rc (r,r\£ n *)tfn*(rV 3 r' = E n k^nk(r). ( ; 


r_. 


1 ) 


with T rc (r, r'. E n k) the non-local and energy dependent exchange- 
correlation self-energy operator. In practice, however, the wavefunctions 
of the two equations are quite close to each other and the effects of replac- 
ing the v TC by can be treated by perturbation theory. The self-energy 
operator can be approximated by Hedin ; s GW approximation [10]. Car- 
rying out this calculation in practice is rather cumbersome and has been 
done for only a few systems [11, 12, 13, 14]. The main finding, however, is 
that the spectrum of E n k and e n k mainly differ by a gap correction which 
is roughly k-independent and energy independent. The latter can actually 
be estimated easily by means of an approximate treatment of the above 
theory proposed by Bechstedt and Del Sole [15]. 

The fact that LDA band structures underestimate the gap is well known 
as the “gap problem". It is actually not a failure of LDA but a question 
of appropriately interpreting the theory. The GW approximation predicts 
the correction to be only weakly dependent on crystal structure. Thus LDA 
band structures e n k can be used to investigate changes in band structure 
with crystal structure. The GW calculations actually also predict a small 
rigid shift of the bands with respect to the average electrostatic potential. 
This is irrelevant for most purposes, because one is usually interested only 
in iuterband differences, but it does have an effect on the band-alignment at 
a heterojunction. The magnitude of this shift appears to depend somewhat 
on the actual choice of parametrization of the LDA exchange correlation 
potential and to that extent is not yet well established theoretically [12]. 
While for wide bands the only important correction appears to be the above 
described gap correction, narrow bands such as core and semi-core levels 
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are subject to additional seif-energy corrections. This is e.g. important for 
Ga 3d states in GaN [16]. 

The method we use to calculate the energy bands, i.e. to solve Eq.(2) ana 
obtain the self-consistent potential, is the linear muffin-tin orbital method 
[IT]. Most of the results reported here were obtained in the so-called atomic 
sphere approximation (ASA) in which Wigner-Seitz spheres replace the ac- 
tual Wigner Seitz ceils and a spherical approximation is made to the charge 
density and potential inside each sphere. This method is in very good agree- 
ment with the so-called full- potential calculations which do not make this 
approximation except for quantities which involve small energy differences 
or in cases where significant distortions of the local spherical symmetry 
occur. For example, cohesive energies, lattice constants, bulk moduli and 
band structures are usually given quite accurately by the ASA. Surface 
energies and workfunctions. however, are not because of the strong lion- 
spherical dipole potential near a surface. In fact, the main problem with 
the ASA is in the spherical approximation to the charge density rather than 
to the potential. The errors induced by using slightly overlapping spherical 
charge densities can to a good extent be corrected by means of the so-called 
Ewald correction [17] in which one corrects for the overlap bv redistribut- 
ing the overlap charge into a homogeneous background charge. We have 
used this approach successfully in our structural energy minimizations for 
the chalcopyrite II-IV-N 2 compounds [18]. For the calculation of optical 
phonon frequencies and elastic constants, which require the calculation of 
very small energy differences in order to stay within the linear response 
regime with respect to the strain or phonon distortion perturbation, we 
need full- potential calculations. The latter were done using the approach of 
Methfessel [19]. 

Our approach for dealing with disorder in alloys and other specifics will 
be described along with the results in the appropriate sections. 

4. Results 

4.1. DIAMOND. BN, AND BEO. 

In Table 2, we present calculated lattice constants, bulk moduli, bond 
strengths (defined as 1/4 of the cohesive energy per zincblende unit ceil with 
respect to free (spin-polarized) atoms), and band gaps for diamond. BN. 
and BeO. We include calculations both for the cubic and hexagonal crys- 
tal structures: i.e. cubic diamond and lonsdaleite diamond and zincblende 
and wurtzite for the other binary compounds. Our results for hexagonal 
diamond and BN are are in good agreement with those of Salehpour and 
Satapathy [21] and Gorczyca and Christensen [22] respectively, but are 
more complete. (Not all quantities given here are provided in the above ref- 
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erences. | We note that the groundstate equilibrium properties are similar in 
t lie hexagonal and cubic forms. The present ASA calculations, however, are 
not sufficiently accurate to determine the relative stability of the wurtzite 
and zincblende structures. In addition, special care has to be taken in us- 
ing equivalent k-point sets and in relaxing the internal degree of freedom 
to establish this small energy difference. That would require full-potential 
calculations. The structure determination was not the primary purpose of 
the present study. We aiso note that typical errors for bulk moduli are 10 
%. The ~10 % difference in bulk modulus between cubic and hexagonal dia- 
mond and c-BN is probably overestimated in part because we did not relax 
the structural degrees of freedom, i.e. the relative position of the two sub- 
lattices and the c/a ratio in hexagonal structure, which are not determined 
bv symmetry. The overestimate of the bulk modulus is consistent with the 
underestimate of the lattice constant and overestimate of the bond energy. 
The band gaps of the hexagonal crystals, however are expected to be quite 
accurate. 


TABLE 2. Properties of diamond. BN. BeO and BeCNj: lattice constants ( a c for cubic, ah 
and c/a for hexagonal) (A), bulk modulus B (GPa), bond energy Eb (eV), and minimum band 
gap Eg (eV). 


cubic hexagonal 




theory 

e.xpt. 



theory 

expt 




Diamond 

a c (A) 

3.58 

3.567 



o* (A) 

2.46 

2.522“ 









c/a 

1.633“ 





B (GPa) 

434 

442 



B (GPa) 

506 





Eb (eV) 

4.2 

3.7 



Eb (eV) 

4.7 





E 9 (eV) 

4.4 

5.5 

r 

- A 

Eg 

3.6 

(4.7)'’ 

r - 

K 

BN 

a c (A) 

3.62 

3.616 



ah (A) 

2.51 

2.56 









c/a 


1.633“ 




B (GPa) 

381 

290-465 



B (GPa) 

425 





E k (eV) 

3.8 

3.3 



Eb feV) 

4.2 





E, (eV) 

4.5 

6.4 

r 

-A' 

Eg (eV) 

5.7 

(7. 6)* 

r - 

A' 

BeO 

o c (A) 

3.80 

3.79 



oh (A) 

2.64 

2.664 









c/a 

1.633° 

1.623 




B (GPa) 

224 




B (GPa) 

248 

224-249 




Eb (eV) 

3.4 

3.0 



Eb (eV) 

3.7 

3.0 




E, (eV) 

6.9 

(9.3) 6 

r 

-.V 

Eg (eV) 

8.3 

10.7 

r - 

r 


“ assuming ideal c/a = ^8/3 

estimate assuming gap correction to be independent of structure 
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The baud gaps depend significantly on the structure. Even the k-pomt 
location of the conduction band minimum changes. In Fig. i we show how- 
some of the conduction-band eigenvalues relative to the valence band max- 
imum change across the series. One may note that the s- Like T[ states in 
both cubic and hexagonal structures come down in the series while the 
r« state of hexagonal and state of cubic move up. Mh. and R\ states 
of the hexagonal structure also move up toward BeO. One can relate the 
k-points of the cubic fee and hexagonal hep Brillouin zones by assuming 
the (lll] c direction to be equivalent to the (0001 ]/, direction and the (110] c 
direction to be equivalent to the ( 1 120]^ direction [23. 24]. The cubic A 
point corresponds to the point 2/3 along M-Lj, and is different from A'/ t . 
We note that in all cases going from cubic to hexagonal the X conduc- 
tion band minimum moves up in energy and the K minimum comes down. 
This is thus a structure related effect. In the case of BeO. however, the 
T minimum is lower in energy than h\ and thus leads to a direct gap 
material. We also note that the gap in hexagonal diamond is lower than 
in cubic while the inverse is true for BN and BeO. The bond energy and 
bulk moduli decrease across the series diamond-BN-BeO while the lattice 
constant increases. Both hexagonal diamond and wurtzitic BN appear to 
be metastable phases. It is not clear whether there is a stability region 
in the phase diagram or whether their occurence is due only to kinetics. 
The latter appears to be the case since hexagonal diamond has sofar only 
been formed under shock compression conditions. BeO is only known in the 
wurtzitic phase. This is consistent with its high ionicitv. The cubic phase 
is significantly higher in energy (0.14 eV/atom according to Chang et al. 
[25]) and is unlikely to be formed. 

Relation to graphitic structures We note that diamond and c-BN have 
competing layered structures graphite and h-BN. which are actually the 
equilibrium phases at ambient pressures. As shown by Fahy et al. [28.29] the 
energy differences between the layered forms and the tetrahedrallv bonded 
forms are small but there is a large energy barrier which impedes the tran- 
sition from one phase to the other. In their models, the transition path 
from rhombohedral graphite to diamond, or h-BN to wurtzite BN corre- 
sponds to a gradual decrease of the distance between the layers which then 
start to buckle. Layered phases of BeO have been shown to be unstable by 
Wentzcovitch et al. [4]. 

Diamond tetrahedral bonding, however, can also be stabilized by sur- 
face hydrogen. Hydrogenation of aromatic molecules yields buckled tetra- 
hedrally bonded rings. Similarly, graphite sheets start to form buckled rings 
when hydrogenated from the edges inwards. This leads to a curling or bend- 
ing of the planes [30] and to the formation of a diamond-like structure on 
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Figure 1. Conduction baud eigenvalues of diamond, BN and BeO versus effective cubic 
lattice constant. Full lines and symmetry notation on the left refer to the hexagonal 
structure, dotted lines and symmetry notation on the right refer to the cubic structure. 


the prism planes of graphite. This is suggested by the open circles in Fig. 2 
which shows the calculated structure at the prism plane interface between 
graphite and diamond. That relaxed structure was calculated 131] using the 
seini-empirical TersofF potential [32]. It was conjectured by Lambrecht et 
al. [31] to form a critical step in the nucieation of diamond. 

Interfaces We have also studied interfaces between some of these materi- 
als: diamond/c-BN [20], diamond/BeO [33] and diamond/Cu [34]. We found 
that the diamond/c-BN offset is of type II. i.e. the vaieuce and conduction 
band of diamond are both above those of c-BN. This also has importance 
for alloys, as will be discussed below. The bonding at these interfaces is 
quite strong because there are no dangling bonds and the valence bands of 
the materials overlap each other in a wide energy range. In contrast, the 
bonding at C'u/diamond interlaces is much weaker because the Cu3d-band 
is narrow and essentially only the Cu4s bands form bonds with the diamond 
dangling bonds. An important quantity is the adhesion energy IF,, which 
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Figure i. Structure of the prism plane interface between graphite and diamond. Three 
layers of diamond almost perfectly match two layers of graphite. Hydrogen atoms (small 
circles) satisfy the dangling bonds of the “non- bonded” diamond layers. The structure 
at the interface (open circles) indicates that the curling of buckled rings forms an initial 
diamond nucleus. 


is given by 


2W a = 7.4 -r 7 S - 7m H) 

where 7 ^, 75 are the surface energies of the materials brought together at 
the interface and 7 ; is the interface energy. For diamond/C'u. the adhesion 
is much weaker than for the interfaces of diamond with c-BN and BeO. 
mainly because the Cu surface energy is so much smaller than that of 
these tetrahedrally bonded materials. Surface energies basically scale with 
cohesive energies because they involve bond breaking. Table 3 gives the 
calculated surface energies and adhesion energies. Details can be found in 
(20. 33. 34]. 

The suitability of c-BN and BeO as substrates for heteroepitaxial dia- 
mond growth was confirmed bv several experimental studies [35. 36. 37]. 
While fully developed epitaxial films could be grown on c-BN [36]. only ori- 
ented growth of particles of size up to about 10 /tin in diameter and with 
slight misorientations of a few degrees has been realized sofar on BeO. This 
is consistent with the better lattice match and adhesion energy to c-BN. 
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TABLE 3. Surface and adhesion energies tin J /m 2 ). 


diamond ideal (110) surface energy 6.5 

c-BN ideal (110) surface energy 5.2 

zincblende BeO ideal (110) surface energy 4.8 
diamond/c-BN adhesion energy 5.4 

diamond/BeO adhesion energy 4.6 

diamond 2x1 (111) surface energy 5.5 

Cu (111) surface energy 2.0 

diamond/Cu adhesion energy 2.3 


4.2. RELATED BE COMPOUNDS: BE 2 C, BE 3 N 2 AND BECN 2 

Besides the elemental IV. Ill- V and II- VI materials of the second row of the 
periodic table, three other compounds with tetrahedral or closely related 
bonding are of interest: BeCN 2 , Be 2 C. and BegN 2 . The main properties 
obtained in our calculations for the first two of these are given in Table 4. 


TABLE 4. Properties of BeCNj and BejC. 



BeCNj 

. Be,C 

Structure 

chalcopyrite 

antifluorite 

Lattice constant a (A) 

3.71 

4.29 (4.342)° 

c/a 

1.96 


11 b 

0.30 


Bulk modulus B (GPa) 

333 

213 

Cohesive energy (eV/atomi 

7.7 

6.0 

Band gap (eV) 

4.2 ( 6.3 > c 

1.15 (2.15) c 


direct T — T 

indirect T — .V 


a expt. value in parentheses 

6 11 is the internal displacement parameter of the chalcopyrite structure (18] 
c With estimated gap correction based on interpolation for related materials. 


Beryllium carbonitride BeC'N 2 is todate a hypothetical material [26, 18), 
which is here assumed to have the chalcopyrite structure in analogy with 
other II-IV-V 2 compounds such as ZnGeP 2 . One may think of it as being 
derived from BN by replacing every other B (group-III) by a Be (group-II) 
or a C (group- IV) atom. This substitution maintains the local stoichiometry 
because each N atom is surrounded by two Be and two C atoms. The 
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• ® O 

Se •: M 

Figure i. C'haJcopyrite structure assumed for the hypothetical material BeCNY- The 
structure represented here is idealized. Bond lengths relax by a relative motion of the N 
sublattice towards C and away from Be. If Be and C atoms are both replaced by B we 
obtain the zincblende stucture of c-BN. 

structure is shown in Fig. 3. The bonding in this material is tetrahedral. 
There are several interesting things about this material. 

- It has a lattice constant intermediate between those of c-BN and BeO 
and may thus form a good bufferiayer for diamond or c-BN growth. 
We note that it is more ionic than c-BN and may thus, as BeO. be 
more likely to avoid layered sp" bonded structures. 

- It has a pseudodirect band gap in the sense that the "directness'' re- 
sults from the fact that the zincblende X point is folded onto T In 
the chaicopyrite structure. A strong mixing of T and A' states, how- 
ever. is expected to occur due to the significant difference between Be 
and C potentials and. in addition, due to the lattice distortions. This 
should lead to a significant dipole matrix element and thus strong 
light-emission. 

- Related materials such as ZnGeP? are promising for non-linear optics 
applications [27] due to the occurence of IV- V bonds. If. as is expected, 
similar non-linear optical advantages of the C-N bonds occur here, they 
would be of particular interest because the band gap is already in the 
UV. 
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Be,C A material with bonding cioseiy related to tetraiiedrai bonding is 
Oe>C. It has the cubic antifluorite structure. The latter consists oi' a fee 
.sublattice of carbon atoms with Be atoms occupying ail tetrahedral inter- 
stitices. In the zincblende structure by contrast, only half of the tetrahedral 
interstices of the fee lattice are filled with a counter ion. We have calculated 
the ground state properties and band structure of this material [38] . The 
main results are given in Table 4. The material is an indirect semiconductor 
with band gap and lattice constant dose to those of 3C-SiC. 


Formation energies Under atmospheric conditions Be^C is unstable to- 
wards formation of Be(OH)^ by reaction with water vapor, or BeO with 
O 9 . This indicates that Be^C is less stable than Be(OH )2 or BeO. Neverthe- 
less. a small amount of Be 2 C formation was detected in growth expriments 
of diamond on BeO [37]. This indicates that in spite of the much smaller 
formation energy of Be 2 C than that of BeO. Be 2 C can form under the very 
reducing conditions of diamond growth. This is significant for the prospects 
of synthesising BeCN 2 . Our calculations predict a much higher energy of 
formation for this material, close to that of BeO. The energies of formation, 
defined with respect to the elements in their standard state and being pos- 
itive in the case of exothermic reactions, i.e. when the compound has lower 
energy than the sum of the energies of the elements, are given in Table 5. 
They were calculated by combining our LDA calculated cohesive energies 
of the respective solids with LDA binding energies of the molecules O 2 and 
N 2 given by Jones and Gunnarsson [39]. In this way, the systematic LDA 
overestimate of bonding is cancelled out. The energy of carbon was taken 
to be that of diamond instead of graphite, the latter being almost equal in 
energy. Zero point motion corrections were not included in our cohesive en- 
ergies because for the present purpose, we only need a qualitative estimate. 
We note that a large part of the stabilization energy of BeCN -2 is due to 
the effective relaxation energy of 1.2 eV/atom. The related material Be 3 N 2 
is included in this table. Its structure is an ordered vacancy structure of 
the antifluorite structure. One may see that per equivalent amount of Be 
atoms, i.e. 1/3 mole of Be 3 N 2 , its energy of 47 kcal is also much smaller 
than that of BeCN 2 . This indicates that under reaction conditions where 
either Be 2 C plus N 2 , or Be 3 N 2 plus graphite, or BeCN 2 could form, the 
last is most likely in a thermodynamic sense. It indicates that when Be 
is exposed to C and N simultaneously, BeCN 2 may form exothermically. 
This calculation, being strictly thermodynamic, of course, does not tell us 
anything about possible reaction barriers. 
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TABLE 5. Exothermic ener- 
gies ot formation lin kcai/molei 
of some Be compounds. 



theory 

expt. 

BeO 

153 

145 

1/2 BejC 

14.5 

14 

1/3 Be, NT 


47 

BeCN? 

156 



4.3. SILICON CARBIDE 

Silicon carbide is unique among the wide-band-gap semiconductors in the 
sense that it exhibits poiytypism. Over a hundred different crystal struc- 
tures called polytypes have been discovered [40]. All are tetraliedraily 
bonded but differ in the stacking of what are the {111} planes in cubic 
SiC- (called -3C-SiC or J-SiC) and the basal planes {0001} in hexagonal or 
rhombohedral SiC, sometimes collectively denoted a- SiC. The polvtypes 
can most easily be described as a periodic structure of twin boundaries in 
cubic SiC. When a twin boundary occurs every 3 (2) layers, we obtain (3) 
((2)) which is the Zhdanov notation [41] for <3H-SiC (4H-SiC). When they 
occur every layer we obtain the wurtzite structure which is purely hexag- 
onal and denoted 2H-SiC or (1). Rhomobohedral polvtypes exhibit more 
complex patterns such as (32) or 15R. 

The relative stability of the various polvtypes has been studied exten- 
sively by Heine et al. [42], Some of the main conclusions of their work is 
that polytypes with bands of three and two layers between twin boundaries 
(e.g. 6H. 4H. 15R) are most stable, followed by 3C SiC. while least stable 
is 2H-SiC. On the basis of their total energy calculations, they proposed a 
model in terms of interlayer interactions. The origin of the poiytypism in 
this context is that the first nearest neighbor layer interactions prefer to 
stack cubically with respect to each other while the second nearest neigh- 
bor interaction favors hexagonal stacking and has almost exactly half the 
magnitude of the first nearest neighbor interaction. The third and higher 
neighbor layer interactions are found to be negligible. Their model explains 
the preference for polvtypes with bands of two and three successive cubi- 
cally stacked layers. It also explains why in surface growth, if complete 3D 
equilibrium is not reached, there is a tendency to form cubic stacking. This 
occurs because for the surface layer there is only one next neighbor layer 
interaction to be considered instead of two in the bulk. Finally,- small dif- 
ferences in the phonon contribution to the free energy were claimed to be 
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the origin oi' different temperature stability regions for 4H anti 6H. Experi- 
mentally, it is found that poiytype selection can be influenced by a number 
of kinetic factors besides thermodynamic equilibrium. Growth on vicinal 
surfaces tends to stabilize the same poiytype as the substrate on which it is 
grown through a ledge growth mechanism. Imperfections on large terraces 
tend to promote independent nucleation and cubic growth. Lower temper- 
atures also tend to lead to cubic growth. All of this is consistent with the 
notion that there must be sufficient mobility of species on the surface to 
move to ledges to continue the substrate poiytype and that otherwise cubic 
nuclei tend to form. In case of cubic growth in the [111] orientation one is of- 
ten faced with the problem of so-called double position boundaries (DPB). 
The latter correspond to the two equivalent orientations of a cubic nucleus 
with respect to the substrate. Stress can induce polytpic transformations 
[43]. Also, ion implantation can influence polytypism after recrvstalLization 
of the damaged area [44]. Finally, doping and growth rates may have an 
influence on which polytvpes form during growth [44]. 

Our own work has focused on the relationship between polytypism and 
properties. Integrated properties such as elastic constants [45] are only 
weakly dependent on poiytype while band structure depends dramatically 
on poiytype. This is because band structure essentially corresponds to the 
formation of electronic standing waves and is thus very sensitive to the 
structure [24]. The changes in band structure also lead to important changes 
in optical reflectivity [46]. 

Fig. 4 shows the dependence of the minimum gap and some specific indi- 
rect gaps on poiytype. The polytypes are here ordered in order of increasing 
hexagonality, i.e. the proportion of layers that are stacked hexagonallv. One 
may see that changes in the position of the conduction band minimum oc- 
cur. While the low hexagonality polytypes have their minimum along the 
hexagonal M — L- axis, including the cubic A’-minimum. 2H-SiC has its 
minimum at K similar to hexagonal diamond and BN. One may also no- 
tice that the specific gaps do not vary monotonously with hexagonality. 
Details on how the band structures are related to each other can be found 
in [24, 23]. 

Fig. 5 shows the calculated UV-reflectivity spectra of some SiC poly- 
types compared to experimental data. We note that a single k-independent 
and poiytype independent gap correction of 1 eV as applied. The latter 
not only provides good agreement for the various peak positions in the re- 
flectivity spectra but also for the minimum gaps given above. This lends 
support to our statements in the section on computational method that the 
gap corrections beyond the LDA are relatively structure independent. 

The changes in band structure between the polytypes leads to different 
effective masses, and hence are important for transport [24]. For example. 
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Figure 4. Band gaps of SiC polytvpes. Squares: experiment, solid circle: thoereticai 
minimum gap: diamonds, upward and downward triangles, specific k-point gaps as indi- 
cated: dashed line, nearly linear T — U gap along U = M — L line of hexagonal Brillouin 
zone. 


it is found that (iH SiC has twelve equivalent minima aiong the M — L axes 
dose to M, which are elongated in the c-direction and with only small bar- 
rier between them. This means that at larger temperature the two minima 
displaced along the c-axis can effectively be thought of as one minimum 
with a very large mass anisotropy. This is consistent with recent measure- 
ments of a large anisotropy in mobilities in OH-SiC [47]. 

We have also studied several interfaces involving SiC. SiC/AIN and 
SiC/BP interfaces [48] will be briefly discussed in the section on alloys be- 
cause of the close relationship among these topics and the new results we 
present here on the corresponding alloys. We also studied SiC/TiC inter- 
faces which are of interest for metallization of SiC and in connection with 
ceramic reinforcement of metals [49]. Finally, we studied so-called inversion 
domain boundaries in SiC [50]. These are planar defects that occur often 
during growth on stepped Si {001}. 
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Figure .5. Reflectivity of SiC poiytvpes. Solid lines experiment, dotted lines theory. 


4.4. BORON COMPOUNDS: BN. BP AND BAS 


Table G presents an overview of boron compounds c-BN. BP aud BAs as 
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TABLE 6. Properties of tetrahedrallv bonded boron compounds: 11 
lattice constant. B bulk modulus. Eh bond enerey, E 0 minimum band 
gap, and conduction-band minimum iocation L,„. 




it (A) 

B (GPa) 

Eh ( e V ) 

E n (eV) 

■^’mm 

BN 

present 

1.62 

381 

3.8 

4.5 

.V 


pp6 

3.606 

367 

3.6 

4.2 

.V 


Expt. “• 6 

3.615 

290-465 

3.3 

6.4 


BP 

present 

4.51 

172 

2.8 

1.35 

A 


pp6 

4.558 

165 

2.9 

1.2 

A 


Expt. 0,6 

4.538 

173-265 

2.6 

2.4 


BAs 

present 

4.79 

130 

2.3 

1.06 

A 


ppc 

4.777 

145 

2.6 

1.2 

A 


Expt.“ c 

4.777 

- 

- 

2.1 d 



Experimental values from Lanaolt-Bornstem rabies (53]. 
b Pseudopoieniial calculation by VVentzcovitch et al. foil 
c Pseudopotential calculation by VVentzcovitch et ai. [52] 
u Our estimate assuming correction same as for BP. Expt. values are uncertain. 


calculated using the ASA-LMTO method. We have included data on BN, 
which given above to facilitate the comparison to BP and BAs. The reason 
why BP in the table below has a slightly underestimated lattice constant 
while the others have sligthly overestimated values is that the calculation 
of [48] for BP did not include the so-called combined correction. For the 
present purposes this is not relevant. Our results are in good agreement 
with the previous calculations by VVentzcovitch et al. [51, 52], although our 
values for the gaps differ somewhat. We note the usual “softening" as we go 
to heavier elements. Due to the very low ionicitv mentioned earlier, the BP 
and BAs compounds have a band structure somewhat similar to Si with 
the minimum along the A axis close to X. Experimental band gaps are 
not well known. We estimate a gap correction of ~1 eV for both BP and 
BAs on the basis of their similarity to SiC and Bechstedt and del Sole ’s 
model [15]. We emphasize that the low ionicitv is expected to reduce polar 
scattering which may give these materials an advantage in mobilities. 

4.5. GROUP-III NITRIDES 

We have performed total energy and band structure calculations for the en- 
tire series of Group-Ill nitrides in both the zincblende and wurtzite struc- 
ture. The two forms are close in energy and can both be stabilized under 
suitable epitaxial growth conditions. This has already been realized for GaN 
and InN but not for AIN. At high pressure the nitrides exhibit a phase tran- 


Q$«Q5NAl PAGE « 
GF POOR QUALITY 



20 


W. R. L. Lamoreciu et ai. 


-it ion 10 the rocksait structure. This is a weil known phenomenon due to 
! lie increase in ionicity under compression and the tendency towards higher 
coordination for ionic compounds. It has been studied theoretically by Go r- 
czyca and Christensen i22! and by Munoz and Kune 1 54] and experimentally 
by Perlin et al. [55]. Volstadt et ai. [561 and Ueno et ai. [57] . Our caicuiateti 
groundstate properties can be found in [23! while full details on the band 
structures are given in [58]. 



Figure o'. Trends in selected band states of Ill-nitrides with atomic number in zincblende 
and wurtzite structure. 

Trends in the band structures. The trends of some eigenvalues with, cation 
in the series of nitrides are shown in Fig. 6. It includes the important eigen- 
values near the gap, i.e. the valence band maximum Tj 5 for ZB or f ]’ and Tg 
for wurtzite. the conduction band state T[ and Jvf for WZ and A'f for ZB 
which “compete” for the minumum. The conduction band T[ states in both 
wurtzite and zincblende are mostly s-iike (actually purely s-like in ZB) and. 
being antibonding states, have important cation-s components. One may 
see that these states show the largest variation with atomic number, and. 
in particular, are monotonicallv decreasing with increasing atomic number. 
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Of course, the above statement requires that we compare the different baud 
structures with respect to the same reference ievei. The natural reference 
level used here is the ASA zero of energy [591 which coresponds to the av- 
erage of the point charge electrostatic potential and is dose to the average 
potential in the interstitial region. This is not sufficiently accurate to de- 
termine band-offsets because that requires consideration of charge transfer. 
But it is adequate for our present qualitative discussion. The reason for 
the downward trend with atomic number is well known and also ocurs for 
example in the series C-Si-Ge-Sn. It is simply due to the fact that s-orbitals 
have a non-zero value at the nucleus and as such feel the —Zjr potential. 
Thus, heavier elements have lower s-ievels with respect to corresponding 
p-levels. Since the state is s-like and the A'f-level in zincblende (or dia- 
mond) is a mixture of s and p-like components and the A'f state in wurtzite 
or lonsdaleite is purely p-like. this explains why diamond. SiC and BN are 
indirect gap materials in either zincblende or wurtzite structure while GaN 
and InN are direct in both. AIN turns out to be indirect in the zincblende 
structure and direct in the wurtzite. Even in WZ. however, the indirect 
and direct gap in AIN are close to each other in energy. Fig. 6 shows that 
the variation of Tf is linear from A1 to In but has a different slope between 
B and Al. This emphasizes again the peculiar role of the second row ele- 
ments. In addition, the Tf 5 state is lower in c-BN than the Tf state. This 
again is a manifestation of the relative energies of p and s-levels and thus 
a result, once again of the absence of a p-like core for 2p wave functions. 
We note that the levels shown are LDA levels which are underestimated for 
the conduction band states as discussed above. Again, however, the LDA is 
adequate for the study of the structurally induced changes in band gap. It 
is noteworthy that for zb-InN the LDA gap is actually negative while the 
real gap is 1.9 eV. The gap in this material is thus almost entirely due to 
the seif-energy correction. 

Ionicitv explains why the bandgap variation is so much stronger in the 
III-N series than in the III- As or III-P series. The reason for this is that the 
high ionicity makes the cation component of the Tf state relatively more 
important. The lower ionicity of GaP compared to GaN. for example, is 
consitent with GaP having an indirect gap. Indeed, the Tf state having 
a smaller Ga (heavy element) component is less sensitive to the nuclear 
potential and lies relatively higher with respect to the Ga p -P, like A'f state. 
By the time we get to GaAs. the anion is also a heavier element and thus it 
also tends to bring the Tf state down with respect to p-like states. Ionicity 
also explains why SiC is indirect in both structures while AIN is only direct 
in wurtzite. The reason is that SiC being less ionic has a relatively smaller 
component of the Si (i.e. "heavier" element) in the conduction band Tf 
state. We note that even for Si. the p-iike Tf 5 state lies below the s-like Ti 
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Elastic constants. We have recently also performed full-potential calcu- 
lations of the zincblende Ill-nitrides including uniaxial strain distortions. 
This allowed us to obtain the full set of cubic elastic constants. From the 
latter, the hexagonal elastic constants can be obtained by a simple tensor 
transformation operation [GO] under the assumption that the bonding is 
similar in both materials. The changes due to the twinning of the tetrahe- 
drons can be incorporated in this transformation [60]. 


TABLE 7. Elastic constants for Ill-uitrides in GPa. (expt. 
values in parentheses i 




AIN 

GaN 

InN 


B 

202 

201 

139 

cubic 

C n 

104 

296 

184 


Cl 2 

152 

154 

116 



199 

206 

177 

hexagonal 

Cu 

398 (345)“ 

396 (296 ) 6 

271 (190)'’ 


C \7 

140 (125) 

144 (130) 

124 (104) 


C, 3 

70 (120) 

64 (158) 

21 (121) 


C33 

468 (395) 

476 (267) 

375 (182) 


C44 

96 (118) 

91 (24) 

46 (99) 


“ Tsubouchi and Mikoshiba [61] 
b Savastenko and Sheieg [62] 


Table 7 gives our calculated elastic constants for cubic and hexagonal 
Ill-nitrides [63]. The deviations from the experimental values are of the 
order of 20 % for AIN — except for the C 13 which is difficult to measure 
— and worse for the other two materials. Particularly disturbing is the 
deviation by almost a factor 2 for the Cu and C33 in GaN and InN. These 
large elastic constants are closely related to the Youngs moduli in the c- 
plane and perpendicular to it and should be more easy to determine than 
the shear moduli. These discrepancies are much worse than typical for LDA 
results and suggest problems with the experimental values. We note that for 
Si [19], Diamond [19], and SiC [45], the elastic constants obtained by this 
method are to within 10 % of the experimental values. Only for AIN are the 
experimental values given here based on sound velocity measurements. In 
the other cases, they are based on analyses of X-ray iinewdiths in powders, 
which is a rather indirect and probably inaccurate method for obtaining 
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elastic constants. The poor quality and small sizes of crystals available at 
t he time of the measurements is probably responsible for these experimental 
problems. 


GaSd and. In4d semicore states A noteworthy aspect of the GaN and InN 
bandstructures is the effects of the Ga3d and In4d semicoreievels. Because 
the lattice constants of these materials are fairly small, the Ga or In atoms 
are close to each other and their semicore d states have significant inter- 
actions. Also, they are close in energy to the deep N2s states. Indirect 
interactions through hybridization with N2s thus further broadens these 
bands. The LDA eigenvalues of Ga3d and In4d actually overlap with the 
N2s bauds. This means that the Ga3d and In4d participate in the bonding. 
It is found necessary to include these d-states as valence bands in order to 
obtain accurate values for the ground state properties (16. 23. 64] . On the 
other hand, in X-ray photoelectron spectroscopy (XPS). the Ga 3d levels 
are found to be situated well below the N2s. This is another manifestation 
of the difference between LDA eigenvalues and the true quasiparticle exci- 
tations of the system. These effects are more pronounced for narrow bands 
and increase with increasing energy separation from the Fermi level. We 
found that the XPS spectrum of Ga3d can be well accounted for by means 
of the so-called ASCF approach [16]. 


Interfaces We have calculated band-offsets among III-N semiconductors 
[65. 66]. The calculations were performed for the zincbiende crystals and for 
the {110} interface. Since the valence-band maxima are similar in nature in 
both strucures, one may expect that valence band offsets would be similar 
for the wurtzite. We can then obtain the conduction band offsets by adding 
the experimental gaps. 


TABLE 8. Valence-band 

AEu = Ev ~ E* (eV) and conduc- 
tion-band A E c = E* — Ef offsets 
with A E’ for zincbiende and A£“ 
for wurtzite. 



t\£v 

A El 

A E? 

AlN/GaN 

0.85 

0.81 

1.78 

GaN/InN 

0.51 

1.07 

1.09 

AlN/lnN 

1.09 

2.15 

3.14 
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Tke valence- baud offset of AlX/GaN measured by XPS fG7i is U.S :r0.3 
e V in good agreement with our results. Interface orientation effects need to 
be studied separately and are expected to be of the order of a few O.i eV . 
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Figure 7. UV-reftectivitv of GaN. The bottom panel shows theory and experiment for 
wurtzite and £ J _ c: the second panel shows the theoretically predicted polarization 
dependence for wurtzite: dashed line E || c, full line E _L c; the top panel shows the 
calculated reflectivity for zincblende. 

Optical properties We have calculated optical response functions for AIN 
and GaN. Our results for GaN are shown in Fig. 7. As for SiC polvtypes. 
we find that good agreement with all major peak positions in experimental 
reflectivity can be achieved with a single constant gap correction or ~1.0 
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■'V. There is. however, a significant discrepancy in absoiuce intensity widen 
increases with increasing energy. A possible explanation Tor this is diffuse 
scattering in the experiment. Our calculations Dredict significant differences 
in spectra between zincblende and wurtzite and significant anisotropy in the 
case of wurtzite. These remain to be confirmed by experiment. For wurtzite 
GaN. our calculated results and the experiments by Rife et ai. [681. shown 
here are in good agreement with those of Olson et al. f69l. A more detailed 
analysis of these spectra will be given elsewhere. 



Figuve S. Calculated tefletivitv of wurtzite AIN: full line. £ J_ c. dashed line. £ H c.An 
LDA gap correction of 1.72 eV was included. 


Our calculated reflectivity for wurtzite AIN is shown in Fig. 8. The 
results are in good agreement with those of Loughiii et al. [70) as far as 
peak positions and general shape of the spectrum is concerned in the energy 
region below ~ 10 eV. As for GaN. our calculations predict a significantly 
stronger reflectivity at high energies than observed in the experiment. 

4.6. ll-IV-Nj MATERIALS 

Ail interesting ciass of wide-hand-gap semiconductors, which has sofar re- 
ceived very little attention, is formed by the II-IV-V 2 materials. These are 
derived from the*corresponding III-V compound by substituting half of the 
III elements by a group II and half by a group IV. This does not completely 
specify the crystal structure. We have performed studies of thes materials 
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assuming the ciialcopyrite structure because tlte latter is a simple super- 
structure of zincblende and is the equilibrium form of II-IV-V? inateriais 
other than nitrides such as ZnGeP?, ZnGeAso, MgSiP 2 . In this structure, 
each group-V atom is tetrahedrally coordinated with two group-II atoms 
and two group-IV atoms. We have studied the compounds with 11= Be or 
Mg and IV= Si and C' and also MgSiPo [IS]. We have already briefly dis- 
cussed BeCN-) along with diamond. c-BN and BeO. Among these nitrides, 
only the silico-nitrides with IV=Si have solar been synthesized. They have 
a crystal structure that is derived from wurtzite in a manner similar to how 
chalcopvrite is derived from zincblende. We have sofar only studied the sim- 
pler but related ciialcopyrite structure. Our main results are summarized 
in Table 9. 


TABLE 9. Properties of II-IV-Vj materials: lattice constant a. bond lengths d. 
bond energy Eh. balk modulus B. band gap E, s . 



a (A) 

djv-v (A) 

d n-v (A) 

E h (eV) 

B (GPa) 

Eg (eV 1 

BeCNj 

5.71 

1.50 

1.71 

3.8 

333 

4.62 (6.3i u 

MgCN 2 

4.11 

1.58 

1.92 

2.5 

>10 

3.69 (5.0) 

BeSiNj 

4.10 

1.T0 

1.76 4 

1.88 

1.76 6 

3.0 

240 

4.47 (5.7) 

MgSiN? 

4.44 

1.77 

1.76 b 

2.08 

2.08 4 

2.8 

195 

4.03 (5.3) 

MgSiPj 

5.64 

2.26 

2.59 

1.9 

75 

1.77 (2.5) 


° Including estimated correction, see Petukhov et al. [18] 
b Experimental values, see Petukhov et al. [18] for original refs. 


Generally speaking (except BeCN 2 ), these materials have higher gaps 
than the corresponding Ill-nitrides, which is probably related to their 
higher ionicity. The folding effects of the chalcopvrite lead to conduction- 
band minima at T. except for MgSiP 2 . Except for BeCN 2 , however, the 
valence-band maxima are displaced from T. Although the gaps are thus 
indirect, except in BeCN- 2 , they are very close to the corresponding direct 
gaps. In Fig. 9 we show the gaps versus lattice constant along with those 
of the related nitrides, diamond, and BeO. The figure also includes a calcu- 
lated gap for AIBN 2 in the (idealized, i.e. unreiaxed) ciialcopyrite structure. 
The lines are only guides for the eve, and should not be interpreted as gaps 
of alloys since the latter may exhibit non-iinear behavior with lattice con- 
stant. One may appreciate from this figure that these materials all have 
gaps similar to those of diamond and c-BN but extend over a lattice con- 
stant range where no binary tetrahedrally bonded compounds occur. They 
thus fill the gap between the lattice constants of the tetrahedrally bonded 
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a (Angstrom) 


Figure 9. Band gaps of II-IV-N? materials and related binary compounds versus effective 
cubic lattice constant a (A). The squares indicate wurtzitic compounds, the circles cubic 
or doseiv related structures. Note that for BeO and AIN the equilibrium phase is wurtzite 
while the zincblende phases indicated by z— are hypothetical. 

materials made from elements of the second row (B. C. N) only and those 
made from mixtures of the elements of the second and third row (SiC. BP. 
AIN). This makes them possibly useful as bufferlayers or substrates for the 
growth of c-BN and diamond. The figure suggests that even higher gaps 
can be expected for the corresponding wurtzitic compounds. 

An interesting question is whether the carbonitrides which sofar have 
not been synthesized can be stabilized. From the above cohesive energies 
and our calculated cohesive energy of SiC and diamond it follows that the 
reactions 

MgSiN 2 + ‘2C — MgCNo + SiC. AH = 1 15kcai/mole (5) 

BeSiN 2 + 2C - BeCN 2 + SiC. AH = -84kcai/mole (6) 

are respectively endothermic and exotherimc. The relative unstability of 
MgCN 2 with respect to MgSiN 2 follows from its larger mismatch in atomic 
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radii. Simiiariy BeCNj has less strain to overcome than BeSiNV The syn- 
thesis of BeC'N 2 , which is the most interesting of this class of materials 
because its properties are closest reiated to those of diamond and because 
it is a direct gap material in contrat to its parent compound c-BN. would 
seem to be synthesizable by the substitution of Si by C' in the existing com- 
pound BeSiN 2 - It remains an open question whether alternative tr-bonded 
crystal structures are possible for these materials. Their higher ionicitv 
which favors higher coordination, however, disfavors such structures. This 
has been argued earlier in this paper to be the reason why such structures 
do not appear for BeO. Since BeCN 2 is intermediate in properties between 
c-BN and BeO. we are presently unable to predict whether it will behave 
like BN and exhibit 7r bonding, or whether it will behave like BeO and not 
show 7T- bonding. 

4.7. ALLOYS 

Several alloy systems among the wide-band-gap semiconductors are of in- 
terest. The Ill-nitrides are known to form continuous alloys. Of these, we 
have sofar only studied the Ga r Alt_ x X system [71]. 

Ga T Al\- x N alloys The baud gaps were found to behave nearly linearly 
with concentration and a cross-over takes place between direct and indirect 
gaps in the zincblende alloys at about 40 % GaN. The relatively low energy 
of formation ( A£/ as 0.02 eV /atom and hence the good miscibility we found 
for this alloy system is expected on the basis of the chemical similarity of 
Al and Ga and the close lattice match. 


.4/ r £i_ x ..Y alloys Alloys with of Ill-nitrides with c-BN suffer from a much 
larger atomic size mismatch. We have not yet performed a systematic study 
of them. Only the 50 % compound with the chalcopyrite structure was 
investigated in connection with the related II-IV-N 2 compounds. Our cal- 
culated cohesive energies for AIBN 2 , c-BN and AIN yield a formation en- 
ergy of A E/ « 2 eV/atom, i.e. two orders of magnitude higher than for 
Al r Gai_ r N. The cohesive energy used here does not include relaxation of 
the bond lengths and is thus an overestimate. Nevertheless, it is clear that 
the Al r Bi_ r N alloy system will have a much lower solubility. The band 
gap we obtained for AIBN 2 is significantly less than the average of those of 
z-AIN and c-BN as can be seen in Fig. 9. The bandgap bowing is conven- 
tionally described by means of a bowing parameter b defined by 


Eg{x) = Eg + A Ej(x - 1/2) - 6x(l - x). 


(7) 
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with E.j — [E.j( 0) + £.,(l)]/2 the average gap and ATT = £-,(1) - E.,(0) 
i he gap difference between the endpoints .-l and B. Using our LDA gaps 
lor c- BN. chalcopyrite AIBNi and z-AIN of respectively 4.73. 3.71 and 3.44 
p V . we had b = 1.5 eV. Including the estimated gap corrections, we obtain 
for these numbers respectively: b.4. 5.0. 4.7. b = 2.2 eV. Bond length relax- 
ations are expected to influence this bowing strongly. The present numbers 
are thus only an indication that bowing is large in this system. A similar 
situation is expected to hold for Bj-Gai.^N alloys. 


ffeterovaient alloys The heterovalent alloys among group IV compounds, 
or. elemental solids, and III- V' compounds present an even more challenging 
problem. We have carried out a study of diamond-c-BN alloys [72] and of 
the alloys of SiC with AIN and BP. All of these are reasonably lattice 
matched alloys systems. We have thus neglected bond-length relaxations 
solar. Details of our calculations for diamond-c-BN were presented in [721. 
Our results on the SiC’-AlN and SiC-BP were reported at the MRS Spring 
Meeting 1993 but have not yet been published, except for a brief mention 
of the SiC- AIN results in [73. 58]. 

First, we discuss the short-range order in these systems and our compu- 
tational approach. Because of the donor and acceptor character of IV- V and 
IV-III bonds respectively, there is a significant energetic advantage in com- 
pensating the two types of cation-anion bonds by placing them in close prox- 
imity and by incorporating them in equal amounts. Charge transfer from 
one to the other will satisfy each bond with two electrons and thus yield a 
gain in band-structure energy [74]. The charging of the bonds that accom- 
panies this, however, leads to a mutual attraction. We thus assume that 
low-energy configurations have nearest neighbor IV- V and IV-III bonds. 
For a similar reason, we exclude cation-cation and anion-anion "wrong" 
bonds completely as these nearest-neighbor interactions would lead to high 
energies. We are then left with a simple pseudobinary alloy problem in 
which we consider only disorder on one of the sublattices and place the 
accompanying anions so as to satisfy the above short-range order charge 
neutrality or compensation rules. In the case of SiC/BP. it is necessary to 
take into account the fact that in BP. B behaves as anion and P as cation. 
Thus we assume Si-P and C-B bonds. The other structures were found to 
have sufficiently high energies that they could be ruled out [48]. 

This does not completely specify the relevant local structures but guides 
the choice of physically plausible structures. We then choose the LI? order- 
ing for 25 % or 75 % and £1 0 ordering for 50 % cation mixtures on the fee 
sublattice of the zincblende structure. These form the basis of the so-called 
Connolly- Williams method [75] in which one considers fee nearest neighbor 
tetrahedrons of the type A n B^^ n with .4 and B the two types of cations and 
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n ranging from 0 to 4. Next, we perform iirst- principles calculations oi the 
total energies and band structures of these mixed ordered compounds. We 
then assume that the properties of the disordered system can be obtained 
as a suitable average over these "basis-states". In principle, one should ob- 
tain the distribution of clusters in the alloy at a given concentration and 
temperature bv the minimization of the free energy. Tlie latter includes the 
energy of formation, the configurational entropy, and eventually the vibra- 
tional free energy and elastic strain energy contributions. This approach is 
known as the cluster variation method. 

For the alloys considered, we find that the energy of formation is nearly 
a parabolic function of the concentration and is high. For 50 % SiC-AIN 
it is about 0.25 eV/atom. This is high and implies that there is only 
a. very limited thermodynamic miscibility. The near paraboiicity implies 
that the Connolly- Williams cluster-expansion model can be mapped onto 
i he ferromagnetic Ising model [72]. This allows as to make an immedi- 
ate and simple estimate of the miscibility temperature, which is given by 
0.816 x 2>\Ef( 1/2)/ kg with kg the Boltzman constant and A.£7(l/2) the 
energy of formation of the 50 % compound. This has a value of ~5000 Iv 
in the present case. Lattice relaxation effects and inclusion of longe-range 
correlations may easily reduce this number by a factor of two. The present 
calculation thus only gives a rough estimate. The experimental phase dia- 
gram [76] of the SiC-AIN alloy system indicates a miscibility temperature 
of ~2000°C. Furthermore, it indicates that alloys preferentially form in the 
wurtzite structure except near the SiC rich end where other polytypes are 
stabilized. In any case, for the present purposes, it suffices to point out that 
these heterovalent alloys are likely to be metastable towards phase sepa- 
ration up to rather high temperatures. Nevertheless, once they have been 
formed by non-equilibrium processes, e.g. from a meit or gas phase mix- 
ture. the phase separation process may be slow because it would require 
massive diffusion. The overall conclusion of this is that thermodynamic 
equilbirum considerations may be rather irrelevant to the actual structure 
of these alloys. Thus, instead of calculating the thermodynamic equilibrium 
distributions of the various types of clusters, we can equally well assume a 
random distribution. 


P n = ^*"(1 -I)*-, (8) 

of the 5 types of basic tetrahedron clusters. We use this procedure to obtain 
a first approximation to the band-gap behavior in the disordered alloys. 
Before presenting those results, we note that inclusion of configurations of 
the type we ruled out from the start would drive up the formation energy 
even higher and reinforce our conclusions. Although such configurations 
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will probably appear in smaii amounts due to random errors during growth, 
they can be considered as defects instead of the average structure. In some 
sense, they are similar to anitisite defects in a simple binary compound. 
Our calculations for the interfaces [481 which considered "wrong bonds" of 
the type Si-Al. C-N. showed that these give rise to defect states in the gap. 

Figs. 10 and 11 show our calculated band structures of the ordered com- 
pounds of SiC-AIN and SiC-BP respectively. The Brillouin zone notation 
follows Bradley and Cracknell [77j. Figs. 12 and 13 show the band gaps of 
the ordered compounds we used as "basis" for the cluster expansion as well 
as the disordered alloy average obtained using a random probability distri- 
bution of the clusters. The bowing coefficients are included in the figure. 
We give both the LDA values and the corrected gaps for the ordered com- 
pounds. We assume here that the correction can be linearly interpolated 
between the known corrections for the end compounds. 

We note that the A'" 6 point where the conduction band-minimum ocurs 
for SiC and AIN is folded onto F in the ordered intermediate compounds. 
For the 50 % and 75 % AIN compounds, however, the valence-band max- 
imum does not ocur at T. This is related to the very flat valence-band 
maximum along the F - X direction in the tetragonal L 1 0 structure which 
correponds to T - M in the simple cubic £12 structure. The valence band 
maximum in both cases occurs at R. Both the minimum gap and the small- 
est direct gap are indicated in Fig. 12. The smallest direct gap, in any case 
is only pseudodirect because it corresponds to a transition form the highest 
valence state at T to the folded X zb minimum. Since the perturbation of 
the parent band structures is strong in these materials, the matrix elements 
for these pseudodirect transitions may be appreciable. The curve for the 
disordered alloy gap is based on the minimum gaps. In SiC-BP. the alloys 
with 50 % or more BP have negative LDA gaps. 

We note that in both cases a very large bowing is present. The band- 
gaps in SiC-BP are predicted to be smaller than the almost equal gaps at 
the end points over almost the entire composition range. As in diamond- 
c-BN. [72] this is related to the type-II offset in this system. Indeed, the 
band-offset at heterojunctions among these materials is shown in Fig. 14 
and indicates that the gap of the overall SiC-BP heterojucntion system is 
smaller than that of the two pure compounds. This is thus proposed to be 
the primary cause of the large bowing. SiC-AIN has a tvpe-I offset, i.e. an 
offset with the larger gap enclosing the smaller gap. Nevertheless, a strong 
bowing is also apparent here. In both cases the origin of this phenomenon 
may be argued to be the large charge transfers taking place between the 
constituent semiconductors. If one thinks of the alloy as a molecular level 
mixture of regions of BP (or AIN) and SiC. it is clear that charge transfer 
processes similar to those found at the hetero junction must occur. To what 
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Figure it. Band structures of ordered compounds of (SiCtalBP). (SiC)(BP) and 
(SiC)(BP)j .The gap is indicated bv shading. 
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Fig tire 12. Band gap versus concentration in zmcblen.de (SiCU-ri AIN It alloys: tri- 
angles. LDA minimum gaps for ordered compounds: pluses. LDA pseudodirect gaps of 
the ordered compounds: squares, corrected gaps for the ordered compounds: solid line, 
corrected gaps for the disordered alloys. 6 is the bowing parameter. 



Figure 13. Band gap versus concentration in zincblende (SiC)i-* (BP)i alloys, triangles. 
LDA minimum gaps of the ordered compounds: squares, corrected gaps for the ordered 
compounds: solid line, corrected gap lor the disordered alloys, b is the bowing parameter. 


extent one can transfer this idea to the homogeneous alloys considered in 
our calculations is not clear. Qualitatively, however, it is clear that similar 
B-Si (N-Si) and P-C (C’-Al) bonds occur and that the balance between 
the charge transfers in these bonds determines the electrostatic dipole off- 
setting the potentials in BP (AIN) units from the ones in SiC units. 

We note that for SiC- AIN alloys of wurtzitic structure, a large bowing 
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Figure 14. Calculated band offsets at SiC/AIN and SiC/BP zincblende heterojunctions. 
The valence band maximum £ v of SiC is arbitrarily chosen as zero. If transitivity is 
assumed the figure also gives an estimate of the AIN/BP offset. 


has indeed been found to exist [78] both for the direct T — T gap and the 
indirect T — A gap. Although our present calculations are for zincblende 
derived alloys, corresponding charge transfer effects are expected for hte 
wurtzite alloys. Thus the experimental results provide indirect confirmation 
of the validity of our results. Work on the wurtzite alloys is in progress. 

5. Conclusions 

In this paper, we have reviewed the electronic structure and related proper- 
ties of a wide range of tetrahedrally bonded semiconductors whose common 
feature is that they include early elements of the periodic table, i.e. Be. B. 
C. N and 0. This survey includes diamond, SiC, the binary compounds 
of formula III-N with the group-III element including B. Al, Ga. In: the 
compounds of formula B-V with the group-V element being N, P. and As: 
and the II-IV-N 2 compounds with IV being Si or C and II being Be and 
Mg. It also includes the II- VI compound of the second row: BeO. The Be 2 C 
and BejNj compounds, which have a bonding type that is closely related 
to tetrahedral bonding, were also included as results for them are relevant 
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for the discussion of the synthesis of BeCN 2 - Bandgaps, lattice constants, 
bond energies or cohesive energies, bulk moduli and, in the case of the ni- 
trides, other elastic constants were presented. The strong covalent bonding 
and the deep N2p levels were shown to be the basic reason for the wide 
band gaps occuring in these materials. The interrelationship between vari- 
ous properties, and, in particular, between crystal structure and electronic 
structure was discussed in detail. Special attention was given to SiC which 
forms in a large number of polytypes. The optical reflectivity of some of the 
most important polytypes of SiC and of GaN and AIN were discussed. They 
provide an important experimental check of the calculated band structures. 
We have revealed the trends with atomic number and ionicity. Related is- 
sues such as the electronic structure of alloys and interfaces were discussed. 
Besides an extensive discussion of previously published material, the paper 
includes new results on SiC-BP and SiC- AIN alloys, new results on hexag- 
onal modifications of diamond, BN and BeO, and on the related material 
Be 2 C. We have analyzed the prospects for synthesis of a new matrial BeCN 2 
on the basis of our calculated energies of formation. 
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Properties of SiC and the Group-III-Nitrides 
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ABSTR ACT: We present a brief overview of the electronic structure studies of SiC 
and the Ill-nitrides performed or in progress in our group. The trends in the 
... band structure across the series and the dependence on polytype are discussed 
as well as the relation between band-offsets at interfaces and band-gap bowing. 


1. INTRODUCTION • • • . 

SiC and the Group-Ill Nitrides form a closely related class of wide-band-gap semiconductors. 

Here we present a brief overview of their electronic structure and related properties. This paper 
is mainly meant as a guide to previous and on-going work carried out in our group. A. full 
account has been dr will be presented elsewhere. Details of the computational approach, can: 
be found in the references. Briefly, we use the linear muffin-tin orbital (LMTO) method (see 
Andersen it al 1987 for a recent description) combined with local density functional theory. 
Most of our work used the atomic-sphere-approximation (ASA). The full-potential version (FP- 
LMTO) was used in some cases as indicated. 

2. HULK ELECTRONIC STRUCTURE . — - - — . ^ 

,..iAYe have calculated the electronic band structure for SiC. BN. AIN. GaN. and InN in both 

the zinc-blende- (zb) and wurtzite ( wz) ..structures (Lambrecht and Segail '1992. 1998a). The 
minimum baud gaps at the experimental lattice constant are summarized in Table 1. 


Table I: Calculated minimum band gaps for wurtzite and zincblende SiC and group- 1 II- nit rides 
with experimental values in parentheses (in eV). The numbers in square brackets are obtained * 
from FP-LMTO. others within ASA. . ' 




• p zb ■ 

Z'G 

-- . - 

Err- 

“ 9 __ j-v rr 

Si( * 

iv. - A', 

1.4 (2.417)“ 

r« - a/, 

2.7 [2.6] 

- 


" ‘ - 


r« - A'z 

2.8 [2.2] (3.30) 6 


BN 

F i s - A i 

4.3 (6.4) c 

r« - h i 

5.7 


-AIN 

I is — A i - 

. £3.4 

ri-r, 

4.6 (6.28)'' 


GaN 

r is - r. 

2.0 (3.5)* 

Ti - r, 

2.4 (3.6) / 


InN 

Eis - r, 

r-0.1 

r 6 - r, 

0.3 (2.05) 3 

- •% \ 


■' Humphreys R G . Bimberg D T aficl Chovke W J 1981, Solid State Coimnun. 39 168 
h Patrick L. Hamilton D R. and Choyke W J 1966. Phvs. Rev. 143 526 
1 Chrenko R M 1974. Solid State Conunun. 14 511 
'' P. B. Perry and R. F. Rutz 1978. Appl. Phvs. Lett. 33739 

f Strife S. Ruan .J. Li Z. Salvador A. Chen H. Smith D J. Choyke W .1. and MorkocJI 
. 1990. /. . Vac. Sci. Technol. D 9 .192. 

'I'**- *'**: ‘ - .irv 1 ’ 

* Dingle R. Sell D D. Stokowski S E". Dean P J. and Zetterstorm R731971. Phvs. Rev. B 3" 
,197 

y Tvagaf V A . Evstigneev A M. Krasiko A N. Andreeva A F. and Malakhov V Ya 1977. 

Sow Phvs. Semicond. 11 1257 

The ASA-LMTO calculated equilibrium lattice volumes, cohesive energies, bulk moduli 
and their pressure derivatives for zb and wz SiC. AIN. GaN and zb-BN were reported by 




ORIGINAL PAGE IS 
OF POOR QUALITY 




I.ambreclii find Segall ( 1 992. 1989 ) and agree wpII with Corczyca and Christensen (1903) and 
ot her calculations referenced there. 

Ue note that, t lie local density approximation (LDA) underestimates the gap by 1-2 eV in 
these systems, in agreement with simple estimates of the self-energy correction introduced by 
I lechst e<| i and Del Sole ( J9XX) and applied to these materials by Lambreclit and Segall ( 1992 ). 
from our calculations ol optical reflectance discussed below, the correction is found to be nearly 
a ((instant shift of the conduction bands and to be independent of polytype. For 211-SiC, we 
note that t he gaps at A and .1 ! are very close within ASA wliile FP-LMTO correctly g : ves the 
minimum gap to lie at l\ . The AT t - state must thus be particularly sensitive to noil-spherical 
coin rilnii ions to the potential. Similar FP corrections are thus expected for the nitrides but. do 
not change the conclusions about the location of the minimum bandgap. 

Alt bought the different polytypes have nearly the same total energies, they have unite 
different band structures, including different gaps and different locations of the conduction- 
band minima. The origin of these changes was explained by Lambrecht and Segall i 1992) by 
((impaling t he hand structures of ZD and YVZ in the same Brillouin zone ( DZ). The conduction 
hand-minima at A’,,.. (.V-/,) move up in the zincblende (wurtzite) structure. This is due to the 
lad t hat A',,.. (A-;,) has lower symmetry in the zincblende (wurtzite) structure. These effects 
change the location of the conduction band minimum in BN. AIN. and SiC but not in CaN and 
InN. I iecau.se in the latter, the minimum always stays at f as a result of t lie deeper potential 
for cation-s orbitals for the higher ^-elements. 


2H - SiC 4H - SiC 



Figure 1: Band structure of 2H and 411 SiC. 

I he band-structure of the intermediate 4H polyfype can be understood in terms of folding 
“fleets of the 211 band structure. The band structures of 211 (using FP) and 4H SiC (using 
ASA) are shown in Fig. 1. The bauds on the surface of the 2II-BZ (L - A — // ) are doubly 
degenerate in 2 II because of the presence of a six-fold screw axis along c. These hands are 
essentially folded onto the M — V — A plane of the 4H BZ. the degeneracy being lifted in The 
lower symmetry 111 structure. The bands now become doubly degerenate on the 411 L — .4 — H 
plane corresponding to the six-fold screw axis with double pitch. 

I lie optical response functions (e^u?). c,(*). and the reflectance) have been calculated using 
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i lie iim/lin-tin orbital basis set for 3C. 211 and 4 II Sit'. The results for 30 and 411 SiC were 
presented in .( Lambrerht rt nl. 1993b) and (Suttrop et al. 1993). Satisfactory agreement 
between theory and experiment is obtained for the major peak positions of 30 and 411 up to 10 
e\ including a constant LDA correction. Preliminary results on the nitrides indicate a similar 
conclusion. 

3. STRAIN EFFEOTS 

A separate investigation using the FP-LMTO method was made of the st rain effects on the 
hnud-st incline of 3C-Si( ' ( Lambrecht et til. 1991a). This study included a calculation of t lie full 
'<•1 i>l elastic constants. The elastic constants tensor is only weakly dependent on the polytype. 
Hence, those for one polytype can be obtained (approximately) from those for another by an 
appropriate Coordinate transformation. We thus obtained good agreement with experimental 
data lor t lie elastic constants of hexagonal or tr-SiC'. no distinction being made here between 
t he various ■hexagohVd jiolytypes. 

Hie st rain effects on the band-structure were reported in the form of deformation potentials 
I Lambrerht ft nl. 1991a). The important deformation potentials, in Pikus and Bir’s notation 
and iii eV. for the hand gap edges are' - . 

’ I'r, l )Jii <1 ‘ tl.\ .V}.; ; ’ • 

1.3 -2.1 -(>.() 2 1.3 S.9 

lo our knowledge the only deformation potential for which there are experimental data 
available is t he acoustic deformation potential Z. for the conduction band minimum at X. It 
was determined to be 22 eV by analyzing the: width of cylcotron resonance lines as a function 
•d temporal lire a n<l assuming. a r x T~ :i / 2 behavior (Kaplan et nl. 1985). The range of the 
data, however, is too narrow to determine this behavior unambiguously. An alternative (and 
even better) fit. to the data, shown in Fig.. 2. is obtained if we assume that in addition to 
• he doFforntalion’ potential scattering, there is a significant impurity scattering, increasing with 
increasing 'temperature* and- approximately described by an inverse lifetime r~ 1 = I.3(T - 
26)l() ln s~ l . Al present, we do not have a well developed theory for the linear dependence of 
i In’ impurity scattering contribution. However, an increasing r~ l with increasing temperature 
would be expected for scattering from ueutraj impurity complexes, a situation which is expected 
lor highly compensated material and likely to be the case in SiC samples available nl that time. 
We note i hat Kaplan et al. ( 1993) reported a new value of 12.6 eV for the deformation potential 
lor samples of higher quality. . ... > i! 
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4. HETEROJUNCTIONS AND ALLOYS 

The electronic structure of alloys was studied using the cluster-expansion approach as de- 
scribed by Lambrecht and Segall (1993c) The valence-band offsets have been determined for 
SiC/AlN and Si('/BP {110} heterojunctions (Lambrecht and Segall 1991b) to be 1.5 and 1.2 
<*V respectively. For SiC/BP I lie ofTset is of type II (staggered). In that case, t he effective gap 
ul the system is smaller than the smaller of the two. Consequently, a large bowing is expected 
for the alloy. Indeed, we find a bowing parameter b = -1.4 eV. with b defined by - ' 

£V.( r) = £>; -f- A£Vy(;r - 1/2) - M 1 - .r). ( 1) 

where E ( , = [Ec;i 0) + £V;-(l)]/2 is the average gap and A£Vy = Ec;( 1) - EVy(O) is the gap 
difference between the endpoints. Even for SiC/AlN. however, we find a large bowing of 3.1 
eV. For SiC/AlN. experimental data were reported for the 2IT poly type (Safaraliev cl ul. 1991). 

In qualitative agreement with our calculations for zincblende. a large bowing of 5.2 eV was 
obtained for the indirect minimum bandgap T — A' as well as for the direct T - F band gap 
(h = 1. 1 e\ .) One may speculate that the large bowing for'SiC/AIN and SiC/BP is related to 
the heterovalent nature of these alloys which leads to important charge re-arrangements and 
hence a large dipole at the heterojunctions and large deviations from a "virtual crystal" picture 
with averaged cations and anions in the alloy. 

We also studied the band-gap bowing in zincblende Al.,Gat-. r N alloys (Albanesi <t ul. 1993). 

In that. case, we determined the cross-over from direct to indirect band gap to take place at 
i = 0.6. The band-gap bowing coefficient in this system was found to be negligibly small 
and slightly upward (i.e. b < 0) if bond-length relaxations are included. For the direct gap 
F ir„. - 1 i,. the nearly linear behavior is in qualitative agreement with the experimental data for 
the wnrlzite F|,. - P| t . gap. which is expected to behave similarly. yrl,. , 
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Center during various stages of this work. E A Albanesi acknowledges support from (A)NICEX J; , 
(Argentina). , . ' " f f 
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